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The properties of the pairing interaction in the shell model framework 
are considered with the aid of the exact numerical solution utilizing the qua- 
' sispin symmetry. We emphasize the features which are out of reach for the 

^ ' usual approximate techniques based on the BCS approach supplemented by 

. the random phase approximation treatment of pair vibrations, especially in 

O ! the region of weak pairing where the BCS+RPA theory fails. For the first 

' time chaotic aspects of the mixing generated by the pairing interaction are 

Q>^ ■ studied. The level repulsion and large information entropy of the eigenstates 

CN , coexist with the absence of thermalization of single-particle motion. The full 

\ spectrum of pair vibration in average displays the spin dependence similar to 

^ ' that for a rigid rotor. 
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I. INTRODUCTION 



O \ Many properties of atomic nuclei are essentially governed by pairing interactions. Here 

^3 



and below we call "pairing" the part of residual interaction which scatters the particles 



' between the pairs of time-conjugate single-particle orbits. The isoscalar pairing, which may 
\ be important for N ^ Z nuclei far from stability, will not be discussed here (at a single 
spherical j-level such T = pairing would require pairs with J = 1), although it also can 
be studied in a similar manner. 

Pairing is often viewed as a simple and the most regular part of the nuclear interaction. 
In a low-lying region of the spectrum it forms a pair condensate which strongly influences 
all nuclear properties [jl],^. According to the standard BCS description borrowed from the 
macroscopic theory of superconductivity, the excitation of the system breaks pairs removing 
them from the interaction domain and blocking the scattering phase space for remaining 
pairs. Another, collective mode of excitation implies the redistribution of the pairs with 
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no breaking. As a result, coherent pair vibrations of the condensate are possible. At some 
excitation energy (or temperature) there occurs a sharp second order phase transition to 
a normal heated Fermi-liquid where the pairing effects are usually neglected, again in the 
spirit of the BCS theory of macroscopic superconductors. However, closer consideration 
of the problem for finite systems such as nuclei reveals that there exist nontrivial effects 
resulting from pairing that persist through the entire shell-model spectrum. In nuclei the 
mesoscopic nature of the system smoothens all transitional phenomena resulting in more 
gradual change of properties. 

Another (usually ignored) aspect of the pairing problem is that, along with the very 
coherent nature of the ground state, pairing brings in a certain degree of incoherent mixing. 
Our goal below is to study the pairing problem in full in order to present a consistent 
picture of the physical effects associated with the pairing interaction. For a broad range of 
energies we study the thermodynamical properties, entropy, thermalization and mixing of 
simple configurations as well as coherent pair excitations, and average moments of inertia. 
A thorough consideration of pairing should precede the detailed study of the interplay of the 
pairing with other components of the residual interaction (the shell model analysis reveals 
unexpected results of this interplay). 

We study pairing in the framework of the realistic spherical nuclear shell model, which 
provides a perfect testing ground for discussing the phenomena of interest. There exist 
a number of approximate methods for curing the well known shortcomings of the BCS ap- 
proach for finite systems. However they cannot be used for our purpose. Pairing fluctuations 
extending far below the BCS phase transition cannot be treated with the usual approximate 
techniques. The spectrum of excited states, including those at high energy, is also out of 
reach for the most of approximations. Having this situation in mind, we consider the prop- 
erties of the exact solution of a nuclear many-body problem with the pairing interaction 
using an algorithm based on the quasispin method ||^. Earlier we have shown |^ that such 
an exact numerical solution is convenient, practically simple and can be combined with the 
new approaches for taking into account other parts of the residual interaction. 
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II. THE MODEL 



We consider a general description of the finite Fermi-system in the framework of the 
spherical shell model. Our space includes a set of j-levels (of capacity Qj = 2j + 1) with 
single-particle energies ej. The method used here for solving the pairing problem exactly 
is based on the quasispin algebras present in the pairing Hamiltonian 

h = J:^3N, + Y.VjrP^Pr- (1) 

j h j' 

Here Pj and Pj are pair creation and annihilation operators of particles on a j-level: 

Pj ~ o ^ (^jmO-jm 5 Pj = (Pj) = ^ ^ '^]m^]m ) (2) 
^ m ^ ra 

where the tilde refers to time-conjugate states, = (— )-'^™ a^- _m, and a^m = —ajm- The 
particle number operator in a given j-level is A^^- = J2m o^jm^^jm- 
Three operators, 



close a quasispin SU(2) algebra; the corresponding classification of the basis is very useful 
in practical calculations. The quasispin formalism used here reduces the pairing problem to 
the problem of coupling of partial quasispins Cj for individual j-levels. We use a coupling 
scheme where the basis states are labeled with the absolute value Cj of the partial quasispin, 
or, alternatively, partial seniority Sj = VLj/2 — ILj^ and particle number Nj for each single- 
particle level j related to the quasispin projection In contrast to the occupancies Nj^ all 
partial seniorities Sj are conserved by the pairing Hamiltonian as well as the full seniority of 
the system s = J2j Sj. The decomposition of the space into blocks labeled by the seniority 
quantum numbers simplifies the problem immensely. Although here we consider only one 
kind of particles, the technique can be extended for full isovector pairing [^-|^. 

As a realistic example we consider the shell-model description of the ^^^Sn isotope. The 
valence neutron pairing is well known to play an important role in tin isotopes P,p!0 . 



The model space consists of five single-neutron orbitals /iii/2, c?3/2, -51/2, 5'7/2, and ^5/2 with 
corresponding single-particle energies -9.76, -8.98, -7.33, -7.66, and -7.57 MeV. There are 
601,080,390 many-body states, among those 272,828 have spin zero and can be partitioned 
into 420 different seniority sets. The dimension of the fully paired set of seniority zero (all 
Sj = 0) is 110. 

In our analysis of the properties of pure pairing we set to zero all non-pairing parts of 
the effective interaction and use the pairing matrix elements {L = and T = 1) from the G- 
matrix derived from the nucleon-nucleon interaction |]ll|, with ^'^^Sn taken as a closed shell 



I I4 . In this derivation the Q-box method includes all non- folded diagrams up to the third 



order in the interaction and sums up the folded diagrams to infinite order [jl3[. In this way our 
input quantities in the Hamiltonian (Q) are completely determined. Calculations with these 
parameters, using direct shell-model diagonalization, as well as the EP+monopole technique 



Q are in a good agreement with experimental data |TJ] . We need to stress however, that our 
main results and the general discussion are insensitive to the precise choice of the pairing 
matrix elements. 

III. LOW-LYING STATES AND PAIR VIBRATIONS 

The pairing interaction of Eq. (||) contains the elementary processes of pair transfer 
from one single-particle orbital to another. Since the pairs are not broken, the interaction is 
only capable of mixing the independent particle configurations within the same seniority set 
{sj}. We refer to the amplitude of such a process as that of a pair vibration. Pair vibrations 
were first proposed by A. Bohr ||15[ addressing the collective nature of low-lying 0^ states. 



A number of other studies followed [|T6[ with a discussion of the collective excitations of 
the pair condensate in the BCS-I-RPA framework. The weakening of the static pairing near 
the subshell closure of ^^^Sn implies strong fiuctuations and mixing of configurations in this 
region that may effect the quality of BCS-based theories in the description of the observed 
states in Sn nuclei. 
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As we mentioned, the total number of the states J'" = O"*" in the valence space of ^^^Sn 
that have no broken pairs (total seniority s = J2j = 0) equals 110. The ground state cor- 
responds to the BCS-like coherent alignment of quasispins in the condensate which produces 
minimum energy and average occupancies Nj which in the BCS limit would be determined 
by the Bogoliubov transformation coefficients, v"^ = Nj, m| = 1 — Nj. The pair scattering 
between the shell model basis configurations is responsible for their complicated mixing. In 
what follows we first analyze the degree of complexity of the resulting eigenstates with the 
aid of information (Shannon) entropy and invariant correlational (von Neumann) entropy 
T7| , p^ j3|l and then try to single out the collective pair vibration modes comparing them with 



the random phase approximation (RPA). 

1. Complexity of the eigenstates 

Information entropy la of the eigenstates |a) with respect to the basis \k) was repeat- 
edly used for characterizing the degree of mixing and delocalization of the eigenfunctions 
expressed in the original "unperturbed" basis |T7|JT9| -p2| . As discussed in [|],^ , the annoying 



basis dependence of this quantity may be advantageous providing an additional knowledge 
about the dynamic interrelation between the reference basis and the eigenbasis of the Hamil- 
tonian. For example, it was shown recently ||2^ that for a shell-model system with random 



two-body interactions, information entropy of the eigenstates expressed in the basis of the 
stationary states generated by the realistic effective interaction does not evolve along the 
spectrum being steadily close to the random matrix theory limit. This means that the wave- 
functions of the system with random interactions are delocalized all over the spectrum of 
the realistic eigenstates. Contrary to that, in the mean-field (shell-model) basis both sets of 
the eigenstates display the regular energy dependence typical for the strong mixing at high 
level density. 

In Fig. 0, the information entropy of all 110 states with s = is shown as a function of 
state energy E^- Information entropy of an eigenstate 
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la- H\Ct?l (4) 

k 

is defined through its decomposition in terms of the basis components |A;), 

a 

Each basis state \k) corresponds to a distribution of particle pairs over single-particle levels 
in such a way that each level has an even number of particles and they are all paired. The 
basis wave functions for 110 different pairwise rearrangements of 16 valence neutrons in 
the ^^^Sn shell-model space are exphcitly built by the action of the operators {P])^^ and 
correctly normalized taking into account the blocking effects. 
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FIG. 1. Information entropy of seniority zero states in ^^^Sn plotted versus their energy. The 
middle panel corresponds to the realistic values of the pairing matrix elements, while for the left 
panel the pairing is enhanced by a factor of 100, and for the right panel the pairing is reduced by 
a factor of 10. 
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Three panels in Fig. |T] show the entropy la as a function of energy Ea for various pairing 
strengths measured with the aid of the overall scaling factor G/Gq in pair transfer matrix 
elements. The results are not symmetric with respect to the centroid of the spectrum since 
the states of s = comprise only a part of all states with spin J = 0; the remaining states 
have seniority 4 and higher. Average information entropy of states in the extreme limit of 
randomness given by the Gaussian Orthogonal Ensemble (GOE) is predicted to be |2ll , p5| ,|3|] 



JGOE = ln(0.482Ar), (6) 

where J\f is the space dimension. In our case of A/" = 110, the GOE limit would give 
-^GOE = 3.97. As seen from Fig. ||, the complexity of states increases with the pairing 
strength and reaches the GOE limit for strong pairing. Already the realistic pairing strength 
provides strong mixing of all states with some reduction towards the ends of the spectrum. 
The asymmetry of la as a function of Ea is related to the asymmetry of the single-particle 
energies in the shell. The highest states mainly correspond to the occupation of the highest 
orbits ("negative temperature"). Because of the relatively large energy gap between these 
levels and the remaining (lowest) unoccupied ones /111/2 and ^3/2 , the pairing interaction is 
less effective in state mixing. 

Being a measure of the delocalization of the eigenstates in an unperturbed basis, infor- 
mation entropy cannot distinguish between a coherent superposition of simple excitations 
characteristic for collective modes and a chaotic combination which emerged from com- 
plicated incoherent dynamics 0. (The distinction would require different tools, such as 
invariant correlational entropy ||18[, see the discussion of Fig. ^ below). Fig. [^demonstrates 



the behavior of information entropy at large pairing strengths. The highest excited state 
approaches the GOE limit very slowly since, being formed as the "antipaired" combination 
at weak pairing, it is shifted up already in the second order of perturbation theory so that 
the further mixing becomes suppressed. The similar effect of the shift down exists for the 
ground state as well but here the collectivity significantly increases the value of entropy. 
The ground state entropy is saturated even above the GOE limit indicated in the figure by a 



thin solid line. For very strong pairing, single-particle energies can be neglected turning the 
problem into the degenerate model case, for Vjj' = const ; the realistic interaction is close to 
this limit. In terms of quasispin, different recoupling schemes can be used to construct the 
basis 1^. The basis states used in EP are labeled both by the partial quasispins and their 
projections appropriate for the weak coupling limit with the distribution of particles close 
to the Fermi step function. The basis of total quasispin and different recoupling quantum 
numbers may serve better in the strong pairing (degenerate) limit. In terms of the particle 
distribution, the strong pairing limit with Vjj' = const results in the equipopulation of the 
sublevels ("infinite temperature"). 
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FIG. 2. Information entropy of selected seniority zero states in Sn as a function of the pairing 
strength. The chosen states are ground state (solid line), 50th state (in order by energy, dotted 
Une) and the highest, 110th, state (dashed line). 

At low pairing strengths (the case of G/Gq = 0.1 in Fig. |I]), the quality of the occupation 
number basis is improved due to the increased role of the mean-field energies. Here the BCS 
theory would show no condensate formation since the spacings between the single-particle 
levels are greater than the pairing strength by an order of magnitude. However, in the exact 
solution, the pair fluctuations still provide a considerable mixing that persists throughout 
the spectrum. 
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IV. RANDOM PHASE APPROXIMATION IN THE QUASISPIN FORMALISM 



2. RPA equations 

Our exact solution obtained by a direct diagonalization in the small space of seniority 
s = can be compared with the standard BCS solution supplemented by the random 
phase approximation (RPA) for the collective modes related to pair vibrations. We start 
with formulating the BCS+RPA approach in the quasispin formalism. For macroscopic 
superconductors (with all partial spins equal to j = 1/2) such an approach was used in the 
classic paper by Anderson ||26|. Since the particle number in this approach is assumed to be 
preserved only on average, one needs to introduce the chemical potential /i. 

The pairing field can be written as a linear functional of one of the quasispin operators 

A,{C} = Y.yn'^f- (7) 

j' 

The operator equation of motion for the partial quasispin components with the pairing 
Hamiltonian are 

[q, H] = -2ze'^q + z[q, A,{£}}]+, (8) 
[q,H] = 2^e'^q-^[q,A,{q}U, (9) 

[q, H] = A,{£|}]+ - A,{/:|}]+, (10) 

where [A, = AB + BA denotes the anticommutator, and we introduced the corrected 
single-particle energies e'j = tj + Vjj/2 — fi. 

The BCS solution corresponds to the static mean field alignment of partial quasispins; 
the spread of the single-particle energies plays the role of the magnetic field along the z- 
axis. The direction of alignment in the xy-plane is chosen arbitrarily (spontaneous symmetry 
breaking); for instance, we take it in such a way that (q) = 0. The static solution can be 
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derived from the diagonal matrix elements of the equations of motion. We denote the static 
amplitudes of £| and £| as Zj and Xj, respectively. Linearizing the equations of motion for 
off-diagonal matrix elements, we obtain normal modes of pair vibrations in the harmonic 
approximation. It is convenient to introduce collective coordinates a and corresponding 
conjugate momenta tt and look for a solution of the form 

£J = + + £^ = yj7r + ..., L] = ^ z^a ^ ... (11) 

(in the harmonic approximation the distinction between the coordinates and momenta is a 
matter of convention). Searching for the excitation in a given nucleus, we need to impose 
the average particle number conservation condition. With the definition (^ of the operator 
the total particle number operator is 

^ = E(§^ + 2/:-)-§ + 2^^- (12) 

In the operator expansion (|1T|) this means that the static and oscillating amplitudes of £| 
should satisfy 

^=? + 2E^.> E%=0. (13) 
^ i j 

The BCS equation for the static part is to be found from the expectation value of the 
right hand side of the ?/-equation (H), 

The normalization of the static amplitudes comes from the condition for the maximum value 
of the total quasispin in the ground state. 



■'f| + 2| = 4— T^f^ + ll^f^l . (15) 



4 V 4 / V 4 
The solution can be written now in the form 
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where the quasiparticle energies are Cj = +'>J^f + ^jl'^J'' ^"^^ S^P equation for = 
Aj{X} is the same as in the BCS theory, 

A. = -EHy^'^- (17) 

The first condition ([13D of particle number conservation determines the chemical potential 
in full agreement with the BCS where the occupation factors are f| = (1/2) [1 — (e^/cj)]. 

Using the harmonic oscillator collective Hamiltonian H = (tt^ + uj'^a'^)/2 with frequency 
uj, we obtain for the RPA amplitudes Xj,yj and zj defined in eq. ([IT]): 

uj^Xj = -2e'^y,+2ZjAj{y}, (18) 
- = - 2^.A,{x} - 2A,z,, (19) 

uj'^Zj = 2Ajyj - 2XjAj{y}. (20) 

The spurious mode of this set, u = 0, gives, either from eq. ([T8|) or from eq. (pO]), the 
solution for yj coinciding with that for the static amplitude Xj, eq. (pIBI). This obviously 
corresponds to the Goldstone rotation of the condensate phase. The spurious mode related 
to the change of the particle number is generated by the operator aJ^jZj, see the second 
condition in Eq. (|1^). However it does not contribute to physical modes since the sum over 
j in Eq. ( PD| ) identically vanishes due to the symmetry of the kernel Vjj'. Looking for the 
non-spurious solutions with u ^ and eliminating the coordinate amplitudes Xj and Zj, we 
come to the equation for the momentum amplitude yj which determines the physical normal 
modes, 

(4e^2 - a;2)% = -e,n,A,{y} + 2Z,A,{t}, (21) 
where the last term contains the same functional A, Eq. (0), taken for 

t,=e'^{2y, + ^A,{y}y (22) 
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This secular equation is valid also for the normal modes generated by pairing with no 
condensate present when we have only the trivial BCS solution, Xj = 0, Aj = 0, the 
vibrational amplitude zj vanishes, and the single-particle energies are equal to \e'j\. With 
the extension to higher orders of the same mapping procedure to the collective variables as 
in Eq. (|n|), one can study the anharmonic effects. However, if the anharmonicity is indeed 
important, as for example in the region of very low frequencies and correspondingly large 
amplitude collective motion, it is simpler (at least in the pure pairing problem) to switch to 
the exact solution. 



3. Example: two-level system with off-diagonal pairing 

As an example of the analytically solved RPA, we consider an interesting particular case 
of two levels with capacities f2i 2 and single-particle energies ei_2 when the pairing interaction 
has only the off-diagonal amplitude V12 = V21 = 9 (in this case the sign of g does not matter). 
With the effective interaction parameters Ai,2 = 5'^i,2/4, the set of the BCS gap equations 
(p!7|) takes the form 

Ai = ^A2, A2 = ^Ai, (23) 
62 ei 

which leads to the exact solution 

X2\2_,/2/2 \2\2_^/2/2 
a2 _ ^1^2 '^l ^2 a2 _ ^1^2 fcj ^2 fr,A\ 

The corresponding quasiparticle energies are given by 



X2_L^'2 X2 I ^/2 

A? + 6^^' '^-^^Ai-f6i^' 



^1 ~ ^1x2 , _/2' ^2 — -^2 X 2 , j2-> K^"^) 

with the useful identity 



6162 = A1A2 (26) 

being valid. The BCS solution collapses, Ai^2 — ^ 0, at the critical coupling strength deter- 
mined by A^Ag = e'i^e2^, or at 
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9 ^ 9c = 



(27) 



The secular equation, see (|2T|) and ( p3D , with the help of the identity (^6]), gives, along with 
the spurious mode uj"^ = 0, the physical root 



The behavior of energies and entropy in the BCS phase transition region are illustrated 
in Fig. ^ for a two-level model and Fig. ^ for the realistic case. First we consider the 
two-level model of the previous section see Fig. |^. In the case of half occupancy and two 
levels of equal capacity, N = Qi = Q2, with energies ±e symmetric with respect to the 
chemical potential /i = , the RPA predicts = 8A^ , see Eq. (pS]). Note that in the case 
of many interacting levels with Vjj' = const the spectrum of normal modes starts with the 
lower value u;^ = 4A^ (threshold of pair breaking). 



^2 ^ 4(^2 ^ ^2 ^ 2e;e^) = 4[(e; + e',)' + Aj + Aj]. 



(28) 



4- Exact solution versus BCS+RPA 
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FIG. 3. Two-level pairing model; AT = fii = = 16, only the off-diagonal pair transfer 
amplitude V12 = V21 = g is taken into account. The upper plot (a) displays the invariant entropy, 
see text, of the ground state for 6g = 0.01, 0.02, and 0.03 averaging intervals. The middle part 
(b) shows the excitation energy of the lowest pair vibration state, thick solid line. The thin dotted 
line approximates this curve with the aid of the RPA built on the normal Fermi ground state on 
the left of the phase transition point at g = 0.25 and on the BCS ground state on the right of the 
critical point. Thick dashed line shows the excitation energy of the lowest state with broken pair 
s = 2. This curve is compared with 2e, thin dash-dot curve, the BCS energy of a two-quasiparticle 
state. The lower panel (c) displays information entropy for the ground state (solid line), second 
pair vibration excited state s = (dotted line), and the lowest s = 2 state (dashed line). 
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The panel (b) in the middle shows the energies of the lowest pair vibration state (thick 
solid line) and the lowest state with one broken pair (thick dashed line) as a function of 
the pairing strength g. The excitation energy of the pair vibration is compared with the 
RPA prediction shown by the thin dotted line. The BCS phase transition occurs, in the 
units corresponding to |ei^2| = ^, 9 = Oc = 0.25, where we see the breakdown of the 
RPA and the vanishing collective frequency. Below this point the RPA is constructed on the 
background of the normal Fermi distribution, while above g = Qc the RPA is built on the 
pairing condensate. The limits g and g ^ oo are described well within the RPA while 
near the phase transition large fluctuations make the BCS+RPA description unreliable. The 



failure of the BCS is well discussed in the literature and studied using simple models p7 
In accordance with this instability, a sharp rise of entropy takes place in the vicinity of 
g = gc showing a dramatic restructuring of the eigenstates when the independent pair basis 
states become strongly mixed. One needs to emphasize here that, even for a very strong 
interaction, the limit of the degenerate seniority model is never reached if the pair transfer 
matrix elements are different for different pairs, in contrast to the model with constant 
pairing when all Vjj' = g =constant. 

Fig. ^(b) addresses the question of the nature of the first excited state. The breaking of 
a pair requires two-quasiparticle excitation energy of 2e that is plotted with a thin dashed 
line, which has in the asymptotic limit of strong pairing the value of twice the BCS gap. In 
realistic nuclear systems unpaired particles may still interact attractively therefore reducing 
this energy even further. The lowering of a two-quasiparticle quadrupole mode is especially 
strong in the case of collective vibrations near the onset of deformation. The lowest state 
of pair vibrations, which has zero spin, is usually much higher in energy. It is possible in 
principle that the state of pair vibrations becomes the lowest excited state. For example 
this can be the case when two single-particle orbitals lie very close in energy making a pair 
transfer energetically more favorable then breaking a pair. This feature is not present in the 
two-level model and the lowest thick dashed line, corresponding to the s = 2 state in the 
plot in Fig. ^(b), is always lower than the thick solid line representing the pair vibrations. 
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FIG. 4. Upper panel: excitation energies of low-lying states in the shell model for ^^^Sn as a 
function of relative pairing strength G/Gq, the lowest pair vibration states with s = 0, solid line, 
and the lowest s = 2 state, dashed line. Lower panel: information entropy of the ground state 
(solid line) and and the most excited paired state s = (dashed line). 
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FIG. 5. Excitation energies of low-lying s = 2 states (dashed lines) are compared with that of 
the lowest pair- vibration state s = shown (solid line) as a function of the relative pairing strength. 
The three panels correspond to ^^^Sn, ^^^Sn, and ^^'^Sn isotopes. 

In the ^^^Sn model, the upper plot of Fig. H, the lowest excited 0+ state is of a pair- 
vibration nature. This pattern persists through the entire range of couplings displayed. The 
same trend is observed in other Sn isotopes, Fig. ||, including the case of ^^^Sn with G/Gq 
above 0.7. Here subshell closure requires higher energy for the s = pair vibration which 
in the limit of weak pairing is just a two-particle-two-hole excitation of the normal Fermi 
ground state. It comes from the fact that there exist close single-particle levels that are 
not fully occupied and thus opened for low energy pair vibrations. In ^^^Sn these are the 
^11/25 d3/2, and Si/2 levels, which are separated by only a few hundred keV. Because of the 
presence of other interactions, this prediction for the first excited state is not fulfilled in 
reality. Still, a number of low-lying 0"^ states is observed in Sn isotopes. The nature of 
these states is complex being influenced by other phenomena such as proton core excitation 
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and shape coexistence [Q. Moreover, the effect of other interactions can be important 
through the admixture of seniority s = 4 states that may have collective structure of a two- 
phonon excitation. In Fig. ^ excitation energies of the two lowest pair vibration states are 
compared with the experimentally observed states that are expected to be mostly of a pair 
vibration nature, see |]2B| and references therein. The agreement is satisfactory considering 
the crudeness of the pure pairing model and the uncertainties in the pair transfer strengths 
and single-particle energies. 
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FIG. 6. Excitation energies of the lowest two pair vibration states are plotted as a function of 
the nuclear mass A for Sn isotopes. Experimentally measured energies of the lowest spin-0 states 
that are expected to have strong s = components, are shown for comparison. 

It has already been discussed in the literature that even though pairing is generally quite 
strong in Sn isotopes there is some weakening in the region of ^^^Sn due to the nearly 2 MeV 
sheU gap between the /iii/2, c?3/2 and Si/2 levels and the rest of the single-particle orbitals 0. 
The proximity of the BCS phase transition results in a rapid growth of information entropy in 
this region, see Fig. ^, lower panel. This also leads to a large excitation energy of the lowest 
pair-vibration state in ^^^Sn at relatively weak pairing. Fig. |^. The large configuration 
mixing in the region of a subshell closure is supported by experimental observations in 
the two-neutron pick-up reaction pD[ and has a direct manifestation in the increase of 



information entropy. 



The randomizing aspect of pair vibrations is somewhat hmited to the region near the 
BCS phase-transition, where many recouphng schemes of quasispin become equiprobable. 
A similar phenomenon for angular momentum is known as geometric chaoticity |^ where 
the many-particle states of good total spin constructed by various recoupling schemes have 
an equivalent mixing (entropy). With no preference to any particular way of coupling, a 
minor change in the interaction may lead to a complete restructuring of the spectrum, a 
phenomenon similar to chaoticity. To emphasize the appearance of chaoticity in the BCS 
phase-transition region the invariant entropy for a two level model is presented in Fig. ^(a). 
The invariant entropy (basis independent) shows the degree of restructuring that happens in 
the eigenstate in response to a change 6g in the interaction parameter g . The invariant (von 
Neumann) entropy of an eigenstate \a) is determined with the aid of the density matrix 

^„ = -Tr{p"lnp°} , (29) 

which is defined as an average over a region of parameters 

Ptk'{9.59) = Yj'^'' CtrCt. (30) 

The invariant entropy in Fig. ^(a) is peaked in the phase transition region, indicating a 
very sensitive (chaotic) dependence of an eigenstate on the interaction strength. Similar 
sharp enhancements of invariant entropy near the phase transitions has been found in the 



interacting boson models |31 . 



V. LOCAL SPECTRAL STATISTICS 

The nearest neighbor spacing distribution P shown in Fig. |^ is very close to that of the 
GOE exhibiting the features of level repulsion. However, the tail of the distribution is not 
Gaussian; it reveals an enhanced fraction of large gaps between the levels. 
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FIG. 7. The nearest neighbor spacing distribution of s = states in Sn shown, fihed his- 
togram; the GOE distribution, open columns. 

The so-called spectral rigidity, or A3 statistics, is the more sensitive probe of spectral 
chaoticity. Fig. | indicates a pseudo-random character of the spectra of the exactly solved 
pairing problem. The spectral rigidity is given as a function of L, the length of the fitted 
interval, whereas averaging over all starting points of such (overlapping) intervals is assumed. 
In Fig. |] the spectral rigidity of paired states in ^^^Sn is compared with the chaotic (GOE) 
distribution and regular Poisson-type behavior. The linear low L behavior of the spectral 
rigidity is Poisson-like, while for larger intervals, starting from around L = 7, the regularities 
of the spectrum flatten the curve as an evidence of the oscillatory behavior associated with 
a regular quasi-harmonic component in the spectrum. In fact, we are dealing here with 
the non-standard situation of a local level repulsion supplemented by the presence of a 
vibrational order. 
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FIG. 8. The spectral rigidity (A3 statistics) of paired states in ^^^Sn compared to that of the 
GOE and Poisson distributions. 



VI. PROPERTIES OF PAIRING THROUGH THE ENTIRE SPECTRUM 

5. Level density 

Even if the mixing is strong the amount of states with a given seniority is very hmited. 
For the lowest seniority this number is small compared to the total number of states in the 
truncated Hilbert space. Globally the level repulsion, Fig. 0, is responsible for the widening 
of the entire spectrum. The density of states in the Hamiltonian systems governed by two- 
body interactions and described by finite matrices is generically close to Gaussian (the same 
is true for the limiting case of no residual interaction ||3^). Pairing is very typical in this 
sense. Fig. ^ shows the density of states in ^^^Sn found from the exact solution of the pairing 
problem. The resulting curve is well fit by a Gaussian. 
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FIG. 9. Density of states in ^^^Sn from the exact solution with the pairing interaction as a 
function of energy. 

The parameters of the Gaussian are the lowest moments of the Hamiltonian which con- 
sists of single-particle and pairing parts, 

H = H,^ + ^H^. (31) 

The centroid of the density distribution is given by the trace of the Hamiltonian, 

H = ^Tv{H), (32) 

where M is the total number of states (or levels for the case of J = 0) in the shell-model 
space under consideration. The width of the level density is the mean square deviation of 
the Hamiltonian, 

a^=lP-H\ (33) 

In general the dependence of the width on the coupling strength G/Gq can contain a linear 
term as well as a quadratic term. As shown in the analysis of the statistical properties of the 
shell-model solutions , the two widths, a combinatorial one for noninteracting particles and 
a dynamical one due to the residual interaction, are added in quadratures. This is possible 



if the two parts of the Hamiltonian (|3l|) are statistically independent, H^p Hp = H^pHp. 
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FIG. 10. The centroid, upper panel, and the width of the density of states, lower panel, in the 
pairing model of ^^^Sn as a function of the relative pairing strength. The empirical curve for the 
width is indistinguishable from the fit of Eq. (|3^ . 

In Fig. the mean value and the width of the level density distribution are shown as 
functions of the relative pairing strength G/Gq. The location of the centroid of the level 
density is linear in the coupling strength, while the width is exactly fit by 



where the single-particle (non-interaction) width, cxsp = 3.01 MeV, and the pairing width, 
Cint = 0.857 MeV, add in quadratures. In the region of the realistic strength the single- 
particle contribution to the width is by far the more dominant. 




(34) 
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FIG. 11. Average occupancies of five single-particle orbitals for all 110 s = states in ^^^Sn as 
a function of energy. 



6. Pairing and thermalization 

A closed mesoscopic system can in principle be analyzed solely in terms of individual 
eigenstates and energy eigenvalues. However, in the region of high level density, a statistical 
description operating with thermodynamical concepts is very appealing. It was shown both 



for complex atoms and nuclei |33,3| that the increase of chaotic mixing with excitation 
energy and level density is in a sense equivalent to the process of thermalization and can be 
treated with the help of a temperature. Moreover, the presence of such mixing is necessary 
in order to make the neighboring eigenstates "look the same" [^] and justify the statistical 
averaging as in a microcanonical ensemble In particular, it turned out that the 
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average occupation numbers of single-particle orbitals can be described by the Fermi-Dirac 
distribution even in the presence of strong interactions, although the single-particle energies 
can be renormalized to become those of quasiparticles. Thus, in the statistical region of 
spectrum we come effectively to the picture of a heated Fermi-liquid. 

However, it is not clear that a specific interaction such as pairing can lead to full ther- 
malization. As shown in Fig. |ll|, the occupancies of single-particle levels derived for the 
eigenstates of the pairing problem vary on average monotonously with energy, therefore 
showing a trend to thermalization. At the same time, there are wild fluctuations which 
would be even further enhanced if states of different seniorities were included. 

To address the thermalization we can define various temperature scales p3|,|3|. The 
thermodynamical "absolute temperature" comes from the density of states 



dE 



- abs 



where p{E) = e 



SpiE) . 



therefore 



TahsiE) = = 



a 



(35) 



(36) 



H-E 

The "single-particle" temperature Ts_p is linked to the termalization of the system in terms 
of the distribution of particles over single-particle levels. It can be determined for each indi- 
vidual eigenstate from the best fit of the occupation numbers to the Fermi-Dirac distribution 



e('=j-A*)/T's-p _|_ X 



(37) 



where the chemical potential fi also can be determined from the two-parameter fit. 
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FIG. 12. Absence of spectral thermalization in ^^^Sn from the comparison of absolute tempera- 
ture, Eq. (35), solid line, with single-particle temperature that follows from the neutron distribution 
over single-particle orbitals, Eq. (37). 



In Fig. Othe comparison of these two definitions of temperature for all states in ^^^Sn 



is shown. The result is very different from what was observed in the shell model taking 
into account all residual interactions where T^bs and T^-p were in a good agreement for 
the majority of eigenstates (a similar thermalization occurs as well when random residual 
interactions are taken). With the pairing interaction only, the dots representing Ts_p of 
individual eigenstates show a low temperature along the entire spectrum whereas the smooth 
hyperbola corresponds to absolute temperature which has a singularity in the middle of the 
spectrum. It is clear that the single-particle thermometer cannot measure correctly the 
temperature of the paired system. 



VII. ROTATION OF A PAIRED SYSTEM 

A pure pairing fermionic condensate has properties similar to those observed in macro- 
scopic superfluid systems. Heating, i.e. the increase of excitation energy, results in a gradual 
destruction of the condensate and an increase in the normal component of the fluid (increase 
of seniority in nuclear physics terminology). Similarly, a higher angular momentum can be 
created only by a pair breaking. The unpaired particles can be coupled to a nonzero angular 
momentum which can be treated as some kind of rotation. When the spins allowed for a 
given seniority are fully aligned, any further increase of angular momentum requires a change 
of configuration and seniority jump, with the corresponding energy increase. This picture is 
qualitatively similar to the phenomenon of quantized vortex formation in superfiuid liquid 
^He. 

In Fig. |T^ the energies of all eigenstates in ^^^Sn within the shell model space (all 
possible seniorities) are marked by points in the energy versus angular momentum plane. 
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The stair-case yrast line, shown by a sohd hne in Fig |T3|, is drawn through the states of 
maximum spin for each seniority. It can be fitted well by a parabola (dashed line), thus 
giving the average moment of inertia / = 32 MeV~^ . This result can be compared with 
the classical moment of inertia of a rigid sphere with uniform density, Jrig = 2MB? /b with 
mass M = 938A MeV and radius R = 1.2A^/^ fm; for ^^^Sn /^g = 35 MeV"^ . It is known 
that in the adiabatic (linear response) approximation a cranked Fermi gas has a rigid body 
moment of inertia while pairing interactions usually reduce this value |l|J^ by about a factor 
of 2. The above results show that the overall effect of pairing on the rotational properties 
is considerably quenched at high excitation energy. The broken pairs and pair vibrations 
return the moment of inertia to the original Fermi- gas value. This phenomenon has a direct 
analog in superfluid liquid ^He, where a large number of vortices emerges at high rotation 



speed forming a lattice structure I^S] which rotates as a whole restoring the rigid body 
moment of inertia and characteristic average velocity field. 

Pairing vibrations are primarily responsible for mixing of states in paired systems and 
therefore they are expected to have a profound effect on the properties of nuclear systems 
such as the high-spin moment of inertia. For calculating the yrast line shown on Fig. as 
a dashed-dotted line, we assume full degeneracy, making all ej equal. This reduces energy 
required for the pair transfer and therefore makes pair vibrations stronger and more chaotic. 
As a result, the moment of inertia was even further increased to a value of 43 MeV^^ . As 
could be expected in the degenerate model, where the perturbative cranking approach is not 
valid, the average angular momentum dependence on energy is vibrational (close to linear), 
and the quadratic contribution related to the inverse moment of inertia is small. 
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FIG. 13. Energies of all many-body states in ^^^Sn versus their angular momentum; the yrast 
line (solid line) is fit by a parabola (dashed line); the dashed-dotted curve corresponds to the yrast 
line in the degenerate case with all single-particle energies being equal. 

VIII. CONCLUSION 

In this paper we discussed the properties of the pairing interaction in a finite quantum 
system described within the shell-model approach. This interaction, being a very important 
part of the general residual interaction has specific features being effective for pairs with 
total spin zero only. The standard approach takes into account the pairing effects with the 
aid of the BCS approximation borrowed from theory of macroscopic superconductivity. This 
approximation fails not only from the viewpoint of particle number nonconservation which 
might be essential for finite systems, but it also does not work for weak pairing and beyond 
the phase transition point. Therefore we base our study on the earlier suggested exact 
numerical solution where the use of the quasispin symmetry provides necessary quantum 
numbers and simplifies the problem immensely. 

The collective excitations in the paired system (pair vibrations) are traditionally de- 
scribed by developing the random phase approximation on the background of the BCS 
ground state. We show that the exact solution differs from that in the BCS+RPA theory. 
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especially in the region of the phase transition. It is seen clearly in the example of the real- 
istic shell model for the ^^^Sn nucleus as well as in the analytically solvable two-level model. 
Analogous results can be obtained for isospin-invariant pairing which will be considered 
elsewhere. 

For the first time we have discussed the chaotic (incoherent) aspects of the pairing inter- 
action. The residual interactions in general mix independent particle configurations creating 
complicated stationary many-body states. These effects influence both local and global sta- 
tistical properties of the system. In the case of a generic residual interaction, we come locally 
to the GOE-like nearest level spacing distribution and enhanced spectral rigidity whereas 
the global behavior reveals thermalization of the system. Typically, the interaction plays 
the role of a heat bath bringing the single-particle occupancies close to the Fermi-Dirac dis- 
tribution even in a strongly interacting system. The pairing interaction is in many respects 
exceptional. Although it induces the local level repulsion and increases information entropy 
of the eigenstates, the spectral rigidity shows a pseudo-oscillatory behavior related to the 
vibrational character of excitations. Both information (Shannon) and invariant (von Neu- 
mann) entropies distinctly reflect the existence of the phase transition. The normal behavior 
of thermodynamic entropy associated with the level density coexists with the absence of full 
thermalization of single-particle motion. The yrast line of the system reveals the trend to 
the rigid-body rotation at high excitation energy where many pairs are effectively broken. 

Pure pairing cannot exist as the only two-body interaction (for a charged system it would 
violate the gauge invariance, see also [0). As shown by the full analysis of the statistical 
properties of the nuclear shell model wave functions 0, the pairing can be considerably 
modified by the presence of other parts of the residual interaction. For example, the fam- 
ilies of states corresponding to various values of the seniority quantum number are almost 
completely mixed although the pairing phase transition is still observed through the behav- 
ior of the pairing correlators for individual low- lying eigenstates 1^^. However, beyond the 
phase transition there exists a long exponentially decreasing tail of dynamic pairing correla- 
tions. The interplay of pairing with other types of residual interaction is an interesting and 
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promising problem for future research. 
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